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Abstract. The resolvent approach is applied to the spectral analysis of the 
heat equation with non decaying potentials. The special case of potentials with 
spectral data obtained by a rational similarity transformation of the spectral data 
of a generic decaying potential is considered. It is shown that these potentials 
describe N solitons superimposed by Backlund transformations to a generic 
background. Dressing operators and Jost solutions are constructed by solving 
a (9-problem explicitly in terms of the corresponding objects associated to the 
original potential. Regularity conditions of the potential in the cases N = 1 and 
N = 2 are investigated in details. The singularities of the resolvent for the case 
N = 1 are studied, opening the way to a correct definition of the spectral data 
for a generically perturbed soliton. 



1. Introduction 

The operator 

C(x,d x ) = -d X2 +d% i -u(x), x = (x 1 ,x 2 ) (1.1) 

which defines the well-known equation of heat conduction, or heat equation for short, 
from the beginning of the seventies [[[], ^| is known to be associated to the Kadomtsev- 
Petviashvili equation in its version called KPII 

K - §uu Xl + u XiXiXi ) XJ = -3u X2X2 . (1.2) 

The spectral theory for the equation of heat conduction with real potential u{x) was 
developed in [|| || , but only the case of potentials rapidly decaying at large distances 
in the x-plane was considered. We are interested in including in the theory potentials 
with one-dimensional asymptotic behaviour and in particular potentials describing N 
solitons on a generic background. In trying to solve the analogous problem for the 
nonstationary Schrodinger operator, associated to the KPI equation, a new general 
approach to the inverse scattering theory was introduced, which was called resolvent 
approach. See @]~|l^l for the nonstationary Schroding er o perator and fl6| fo r the 
Klein-Gordon operator. Some results for the operator ( |l.l| ) were given in ||17||. For 
the specific case of N solitons on a background for the KPI equation see especially 

Here we apply the same approach to the heat equation and construct a potential 
u'(x) describing N solitons superimposed to a generic background potential u(x). 
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Superimposition is performed by means of a rational similarity transformation of 
the spectral data of u(x). This procedure supplies us not only with the Backhand 
transformation of the potential u(x), but also with the corresponding Darboux 
transformations of the Jost solutions and with the spectral theory for the transformed 
potential u'(x). All related mathematical entities such as the extended resolvent M', 
the dressing and dual dressing operators v' and w' and the Jost and dual Jost solutions 
and 'J'' corresponding to u'(x) are given explicitly in terms of the same objects 
associated to the background potential u(x). In |l9|, some preliminary results were 
presented for N = 1 and 2 by using recursively the binary Darboux transformations. 

We show that the main mathematical object of the theory, i.e. the extended 
resolvent M', is given as a sum of two terms. The first one is obtained by dressing 
with the operators v' and uj' the resolvent M of the bare heat operator £ (x, d x ) = 
—d X2 + d xi , while the second one, m', takes into account the discrete part of the 
spectrum. The dressing operators v' and lo' are constructed by solving a 9-problem 
for the transformed spectral data. 

The theory with respect to the nonstationary Schrodinger equation is in some 
respects simpler and in some other respects unexpectedly more difficult. We give 
the explicit expression for and u' for any N, but the reality and regularity 

conditions for the potential are rather involved and we examine in details only the 
cases N = 1,2. We study the singularities of the resolvent in the case N = 1. This 



resolvent can be used for investigation of the spectral theory of operator (|l . l| ) with the 
potential being a perturbation of the potential u'(x) obtained by adding to it a 'small' 
function u 2 ( x ) rapidly decaying on the x-plane. It can be shown that under such 
a perturbation of the potential the Jost solutions get singularities more complicated 
than poles, but on the other side they have no additional cuts on the complex plane of 
spectral parameter, in contrast with the nonstationary Schrodinger case. This means, 
however, that also in the case of the perturbed heat equation the standard definition 
of spectral data for a generic non decaying potential must be modified. The solution 
of this problem is deferred to a future work. 

2. Direct and inverse problems in the case of rapidly decaying potentials 

In the framework of the resolvent approach we work in the space S' of tempered 
distributions A{x,x']q) of the six real variables x, x' ', q <E R 2 . It is convenient to 
consider q as the imaginary part of a two-dimensional complex variable q = qjj+iqQ = 
(li j Q2 ) € C 2 and to introduce the "shifted" Fourier transform 

A(p;q) = Jdxjdx' e^ + ^-^ x ' A(x,x';^) (2.1) 

where p £ R 2 , px — p\X\ + p-ix-i and qsRX = qisjfiEi + q2SR a; 2- We consider the 
distributions A(x,x';q) and A(p] q) as kernels in two different representations, the 
^-representation and the (p, q)-representation, respectively, of the operator A(q) (A 
for short). The composition law in the x- representation is defined in the standard 
way, that is 



{AB)(x, x'; q) = Jdx" A(x, x"; q) B(x",x'; q). (2.2) 
Since the kernels are distributions this composition is neither necessarily defined for 



all pairs of operators nor associative. In terms of the (p, q)-representation (2.1) this 
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composition law is given by a sort of a "shifted" convolution, 

(AB)(p;q) = fdp'A(p- p';d+ p')B{p';q). (2.3) 

On the space of these operators we define the conjugation A* and the shift A^ for the 
complex parameter s € C 2 , that in the x- and (p, q)-representations read, respectively, 



A*(x,x';q)=A(x,x';q), A* (p; q) = A(-p; -q) (2.4) 
where bar denotes complex conjugation, and 

A (s) feq) = A(p;q + s), A^(x,x / ;q) = e is ^ x ">A(x,x';q + S ^). (2.5) 

For any operator A, thanks to the fact that its kernel belongs to the space <S'(R 6 ), we 
can consider its dj -differentiation (in the sense of distributions) with respect to the 
complex variables q, (j = 1,2): 

(B j A){p;c l ) = ^^-, (5 j A)(x ) x';q)= 1 -(x j -x' j + -£-)A(x,x';q), (2.6) 



where the formula for the kernel of djA in the x- representation is obtained by (2.1). 

The set of differential operators is embedded in the introduced space of operators 
by means of the following extension procedure. Any given differential operator 
A(x,d x ) with kernel 

A(x,x') = A{x,d x )8{x-x'), (2.7) 

5{x — x 1 ) = 8{x\ — x' 1 )5(x2 — x' 2 ) being the two-dimensional <5-function, is replaced by 
the operator A(q) with kernel 

A(x, x';q) = e^-^ A(x, x') = A(x, d x + q)5{x - x'). (2.8) 
We refer to the operator A(q) constructed in this way as the "extended" version of the 



differential operator A. It is easy to see that in terms of the (p, q)-representation (2.1) 
the dependence on the q- variables of the kernels of these extensions of the differential 
operators is polynomial. In particular, let Dj denote the extension of the differential 
operator d x (j = 1,2), i.e. according to ( |2.8| ) 

D j (x,x';q) = (d x .+q j )5(x-x'), j = l,2, (2.9) 

then Dj in the (p, q)-representation takes the form 

r> i (p;q) = -iq J -(5(p). (2.10) 

An operator A can have an inverse in terms of the composition law ( |2.2| ), ( |2.3[ ), 
say AA~ X = I (in general left and right inverse can be different), where / is the unity 
operator, 

I(x,x';q) = 5(x-x'), 7(p;q) = 5(p). (2.11) 

In order to make the inversion A -1 of an (extended) differential operator A uniquely 
defined we impose the condition that the product (A^ 1 )^ A^ 1 exists and is a bounded 
function of s in a neighborhood of s = 0. Let us consider as an example the operator 
D\ — a, where, being a a complex constant, we write for shortness a instead of al, 
according to a general notation we use in the following. For its inverse operator we 
get 

(D 1 - a)- x (x, x'; q) = sgnfe - a^e^^* 'M(Qi ~ " <))^2 ~ 4), (2-12) 
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that is just the standard resolvent of the operator d xi . In terms of the (p, q)- 
representation this inverse operator is given by 

(£>i - a)" 1 (p; q) = i(qi - ia)" 1 ^). (2.13) 

The simplicity of this formula makes clear the usefulness of the (p, q)-representation. 
Notice also that in this representation the necessity of requiring the boundedness 
condition introduced above for the inverse is especially evident since it is just this 
condition that excludes the presence of additional terms of the type <5(q! — ia) in 

According to this general construction the extended operator L(q) corresponding 
to the operator (1.1) is given by 

L = L -U, (2.14) 

where Lq = —D 2 + D\, i.e. the extension of C(x, d x ) in the case of zero potential with 
kernels 



L Q (x,x';q) = I - (d X2 + q 2 ) + (d Xl + <7i) 2 J 5(x-x'), L (p;q) = (iq 2 - qf)<J(p)(2.15) 
and U can be called the potential operator since has kernels 

U(x,x';q)=u(x)6(x-x'), U(p; q) = v(p), (2.16) 
where v(p) = (27r) -2 Jdx e lpx u(x) is the Fourier tran sform of the potential u(x). Below 



we always suppose that u(x) is real, which by ( |2.4j ) is equivalent to 

L* = L. (2.17) 

The main object of our approach is the (extended) resolvent M (q) of the operator 
L(q), which is defined as the inverse of the operator L, that is 

LM = ML = I. (2.18) 

In order to ensure the uniqueness of M we require that the product M^M is a 
bounded function in the neighborhood of s = 0. Then, in particular, the resolvent Mo 
of the bare operator Lq has in the (p, q)-representation kernel 

M ( P;q ) = -^L I . (2.19) 

lc l2 Vll 

In the ^-representation we get 

M (x, x'; q) = 6 ^ J da [9(qf -q 2 - a 2 ) - 6{x 2 - x' 2 )] e -^+^)(— '). (2 .20) 

The resolvent of L can also be defined as the solution of the integral equations 

M = M + M UM, M = M + MUM . (2.21) 

Under a small norm assumption for the potential we expect that the solution M is 
unique (the same for both integral equations) and that it satisfies the boundedness 
condition at s = for M (s) A/. 

The resolvent is directly connected with the Green's functions of the operator 



(1.1). Indeed, since by definition (|2.8[ ) the product e q ( x x ^L(x,x';q) is n othing but 
the kernel L(x,x') of the original operator C(x,d x ) in (1.1), we have from ( 2.18 ) 

£{x,d x )(e q{x - x ' ] M(x,x' ';<?)) = £. d (x',d x ,)(e^ x - x '^M(x,x';q)) = S(x - x'), (2.22) 

where C d is the operator dual to C. Notice that while the product e q ( x ~ x 'L(x,x';q) 
is ^-independent the same combination for the resolvent (see the simplest example in 
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( |2.20 )) essentially depends on q. This means that the resolvent can be considered as 
a two-parametric (q € R 2 ) family of Green's functions of the operator C. 

Thanks to ( gj| ) and the kernel M(p;q), like M (p;q), is singular for 

q = £{qi) and q + p = £(qi + pi), where we introduced the two component vector 

t(a) = (a,-ia 2 ), (2.23) 

such that L (p;£(qi)) = 0. A special role in the theory is played by the operators v 
and lo, whose kernels in the (p, q)-representation are given as values of (MLo)(p;q) 
and (L M)(p;q), respectively, along these curves: 

^;q) = (ML )feq) cj(p; q) = (L M)(p; q) . (2.24) 

q=*(qi) q=^(qi+pi)-p 

It is clear by construction that both kernels v(p; q) and u(p; q) are independent of q2 

and it is easy to see that both of them tend to 8{p) as qi — > oo. In the x- representation 

these kernels are given by means of the inversion of ( ]2.l| ) as 

(2.25) 
(2.26) 



Kzy ; q 9 ) = ^ 4) / d qm/ d ^ ei£!R(qi)(3: '" 3: " ) ( ML o)(^^ , V a (qi)) ; 



2vr 



c^,^) = ^ X>2) jd^Jdx"^^"-^ (L M) (x" , x' ; f g, (q x ) ) . 

Operators f and cj obey the conjugation properties 

i/*=i/, cj*=w, (2.27) 



which arc equivalent to (2.17), and they are mutually inverse, 

ujv = I, vuj = I. (2.28) 

The most essential property of these operators is that they dress M and L by 
the following formulae 

M = i/MqU, L = vLqlu, (2.29) 

that thanks to ( 2.28 ) give 

Lv = vL , ujL = L lj, Mv = vM q , ujM = M u). (2.30) 

Thus it is natural to call v and w dressing operators, or more specifically the Jost 
dressing operators, since the Jost solutions can be directly related to them by means 
of the following construction. Let us introduce 

$(s, k) = e-^'^it, k), *(.t, k) = c^^x, k), (2.31) 

where 

x(x,q 1 )=Jdpe' i P^(p;< l ), £(x, qj = /dpe" 1 ^; q - p) (2.32) 

and where we named k the spectral parame ter q i in order to meet the traditional 
notation. Then the first pair of equalities in ( [2.30 ) takes the form 

C(x, d x )<f>(x, k) = 0, C d (x, d x )$(x, k) = (2.33) 

so that $ and \& solve the heat equation and its dual. In order to prove that they are 
the Jost solutions we first note that for \ an d £ in ( [2.32 ), thanks to ( |2.l[ ) and (2.25), 
( |2.26D , we get 



=yd,'(M £o )«.«.(x,«V a (k)) = /d*'„«(*,*<;0), (2.34, 
£(x',k) =Jdx(L M)^W\x,x';ez(k)) ^Jdxu^Xx^x'-Q). (2.35) 
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The first equalities here can be used in ( |2.24| ) and then we have by (2.31) and (2.32) 
that 



$0, k) = Jdx' (£$(x', d x ,)G(x, x',k))e 
*(x',k) =[dx (C (x,d x )G(x,x',k))e ie ^ x , 



i£(k)x' 



where the Green's function of the Jost solution 
G(x,x',k) = e q ^- x ' s >M(x,x';q) 

appeared. Being a special case of the object considered in ( |2.22| ) it obeys 
C(x, d x )G(x, x', k) = C d (x', d x ,)G(x, x', k) = S(x - x'). 



(2.36) 
(2.37) 

(2.38) 

(2.39) 



From the integral equations (2.21) we have ML = I + M UML and L M 
I + L MUM that thanks to ( [2.36 ) and Q2.37 ) give the integral equations 



k) = e _u W* + J dx' G (x - x', k)u(x')$(x', k), 
*(x, k) = e i£(k) * + fdx' G (x' - x, k)u(x')V(x', k), 



(2.40) 
(2.41) 



where the Green's function Go(x — x',k) is given by means of (2.38) in terms of the 
bare resolvent Mo- It is easy to see that thanks to ( 2.20| ) 

G (x, k) = ±-Jda [0(kl - a 2 ) - 6(x 2 )] e -^"+^>>, (2.42) 

which is the Green's function of the standard integral equations defining the Jost 
solutions. Self-conjugation property of the operators v and co (2.27) can be 
reformulated in terms of the Jost solutions &(x,k) and ^(x,k) in the following way 

$(x,k) = $(x,-k), *(x,k) = ®(x,-k) (2.43) 



as well as relations ( 2.28 ) that become the scalar product and completeness relation 
for the Jost solutions 



JdxiV(x,k+p)$(x,k) =2ir6(p), peM, 
J 6k^(x',k)^(x,k) = 2nS{x 1 - x[). 



(2.44) 
(2.45) 



Correspondingly, we call the first and the second relation in ( [2.2S ) the scalar product 
and the completeness relation for the dressing operators. 

The dressing operators themselves can be reconstructed by means of the Jost 
solutions, 



,( Xi x '-q) = ^2 — ^) e -9i(»i-»l) f dk e U(fc+i9i)«'$( a . ) k + iqj, 
2n J 

o(x,x';q) = 6{X2 ~ X ' 2 K -^-^ 
2n 



Idk e-^+^^ix'^ + iqJ. 
They can be expanded into the formal series 

oc oc 

V = I + Y^V-n{2D i r n , c = / + ^(2 J D 1 )-"^„, (2 



(2.46) 
(2.47) 



n=l 
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that in the (p, q)-representation have the meaning of asymptotic expansion in powers 
of 1/qi for v(p\q) and of l/(qi + Pi) for w(p;q), i.e. the asymptotic expansions of 
the Jost solutions with respect to the spectral parameter at infinity. For the first 
coefficients in the expansion we get 

!/_i+w_i = 0, [D 1 ,v_ 1 ] = -[D 1 ,u_ 1 ]=U, (2.49) 

that reconstruct the potential in the standard way, that is u(x) = d xi u^i{x). 
In order to formulate the Inverse problem we introduce the operator 

(2.50) 

Then taking into account the 



L ML 



that can be considered as a truncated resolvent. 



definition of the derivative in (2.5) it is easy to show that the (^-equations for the 
dressing operators are given in the form 

d x v = vR, B x uj = -Roj, (2.51) 
where the operator R is given by 

R(p;cO = r(c Ll )6(p + 2e Ui (q 1 )), (2.52) 
with the spectral data r(qi) defined by means of the following reduction of p: 

r(q x ) = 2 7 rsgn(q m )/j(p;^(q 1 )) 



D=-2f H (q 1 ) 

Let us mention that by construction the operator R obeys the conditions 
[L ,R} = 0, [M ,R}=0, 
R* = -R. 



(2.53) 

(2.54) 
(2.55) 

(2.56) 

where r(k) = — r(— k) (cf. ( |2 . 5 5| ) ) . In orde r to get th e repr esentation of the s pectral 
data in terms of the Jost solutions we use ([yj) and (|2.18[) and rewrite (|2.50|) in the 
two equiva lent f or ms p = U MLp and p = LqMU . Then for the value of p in ( 2.50 ) 
we get by (|2.24|) , ( |2.3l| ) and (H||) that 



In terms of the Jost solutions equations (2.51) take the standard form 
d k $(x, k) = $(x, -k)r(k), d k *(x, k) = -k)r(-k), 



r(k) 



dxe- ie( - k)x u{x)»f{x,-k) 



(2.57) 



2vr J v ' v ' ' 2ir 

Finally, the normalization conditions that complete the formulation of the inverse 
problem for the Jost solution or its dual are given by (2.31) in the standard forms 
limk^oo x( x j k) = 1 or, correspondingly, lim^oo k) = 1 as follows from (2.32) 



and (2.48) 



3. Similarity transformations of the Spectral Data 

3.1. Rational similarity transformations 

We study transformations of the above introduced objects generated by similarity 
transformations of the spectral data of the form 

R' = WRW' 1 (3.1) 

where the operator W is the following rational function of the operator D\: 



N 



Di - a, 



(3.2) 
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with dj and bj some complex parameters. Such W guarantees that R' is of the form 
( |2.52D with 



'W = IT r ■ J rto ( 3 - 3 ) 




substituted for r. It is easy to see that such R' obeys properties ( |2.54 ), i.e. it commutes 
with Lq and Mq, and in order for R' to obey property (2.55) we have to impose the 
following conditions on the parameters: 

a 3=<MV b i = K(j) ( 3 - 4 ) 
where n a and 7T& are some permutations of the indices. Notice that in the simplest 
situation where N — 1, i.e. W — (£>i — a\)/{Di — foi), the parameters must be 
real, a\ — a\ 1 b\ = b\. Like in the nonstationary Schrodinger case the potential 
u' corresponding to this simplest situation can be obtained from the potential u by 
means of a binary Darboux transformation suggested in [ p2| . However, in contrast 
with the nonstationary Schrodinger equ ation, the generic case, with W given by (3.2) 
and Oj-j bj subjected to conditions ( |3.4|) but not necessarily real, cannot be obtained 
by applying recursively binary Darboux transformations or one needs to admit non 
real intermediate potentials in the iterative procedure. Here we shall not use the 
techniques of the Darboux transformations but rather derive these transformations 
solving the inverse problem given by the spectral data R' explicitly in terms of the 
objects corresponding to the generic original spectral data R. 

Let us mention that the new spectral data R' have additional discontinuities if 
compared with R. Indeed, r'(k) given in ( |3.3| ) is undetermined at points k = ia,j 
and k = ibj while r(k) as given in ( 2.57 ) is a well defined function for kjf ^ 0. 
Correspondingly, these additional discontinuities lead to new properties for the new 
dressing operators v' and lo' . Indeed, we show below the solvability of the (^-equations 

N N 

d x v' = v'R' + iTrJ2 L 'a j S ( D l- a j)> 8^' = -R'uj' + m S ^5(D 1 -b j )uj' bo , (3.5) 



that differ from equations (2.51) by additional ^-terms at D\ — aj and at Di = bj. In 
other words we show that v' and oj' can have poles on the complex plane with residua 
v' a . and u>' b , . Different choices of the form of these additional terms correspond to some 
renormalizations of the dressing operators and, consequently, of the Jost solutions^]. 
One could also consider the case in which poles are absent, but we will show that, 
then, the solution does not contain solitons and decays at large distances. 

In order to close the formulation of the inverse problem for v' and lo' we have 
to impose normalization conditions on v' '(p; q) and to'(p; q) at some value of qi (both 
of them, as we know, must be independent on the variable q2). Taking into account 
that v and lo are normalized by the asymptotic condition at infinity, it is natural to 
choose the same point for the v 1 and lo' . Then it is easy to see that, without loss of 
generality (omitting the uninteresting case of a potential u(x) shifted by a function of 
X2 only), we can fix that 

v'{p\ q) -> 6{p), Lo'(p; q) -> 6(p), qi -> oo. (3.6) 

* In [Q, where the direct problem was examined, in order to define Jost solutions via an integral 
equation invariant in form, i.e. not depending on the parameters aj and bj, a different normalization 
was necessary. This is another intriguing characteristic of the heat equation in comparison with the 
nonstationary Schrodinger equation. 
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In this discussion we used that di(D 1 — a) 1 — iir5(D 1 — a). The kernels of this 
5-operator in the x- and (p, q)-representations are given by 

(5(D 1 -a))(p;q)=6(q 1 -ia)6( P ), 

e ia a (x 1 -x' 1 ) 

(S(D 1 - a))(x,x';q) = — d(q 1 - (h R .)5(x 2 - x' 2 ), (3.7) 

where we used the definition 8{z) = 5(z^)S(zq) for any z € C. Note that in the 
(p, q)-representation z/(p;q) and a/(p;q) have poles, respectively, at qi = iaj and 
qi = ibj — pi and we have 

^LO;q) = -i res v'ip;^), ^ (p;q) = -i res w'feqj. (3.8) 

J qi=iaj J q 1= ib j -p 1 

Correspondingly, the Jost solutions $'(cc,k) and ^/'(a;,k) have poles, respectively, at 
k = iaj and k = ibj and in the x-representation we have 

u' a .(x,x';q) = S(x-x')e^ ia ^' a .(x), (3.9) 

J h . (x, x'; q) = S(x - x')e-^>'^ bj (x') (3.10) 

where 

& (x) = -i res $'(x,k), ^.(x) = -i res *'(x,k). (3.11) 

J k— la^ J k— ibj 

Note that in the formulation of the inverse problem the values of the residua at 
poles v' a . and u>' b are left free. We show below that in order to close the construction 
of v 1 and lu' it is necessary to impose that they dress the Lq operat or, th at is, that 
they obey the differential equations like the first pair of equations in (2.30) 

LV = u'L , uo'L' = L u', (3.12) 

where 

L' = L Q — U' (3.13) 

with some new potential u' . More precisely, we show in Section 3.2 that if v' and u/ 
are dressing operators they satisfy the orthogonality relation wV = / and, then, in 
Section 3.3, that this additional requirement imposed to v' and oj' obtained from the 
solution of the d problem (|T^) is sufficient in order to guarantee that they satisfy the 
dressing equations (3.12). If poles are absent v' and uj' are uniquely determined by 
( |3.5| ) and their asymptotic behaviour at large q and it is not necessary to impose any 
additional requirement in order that they satisfy equations (B.12). 

Let us note that the new dressing operators, since they obey ( [3.12] ), have the 



— n : ^ — n 



same (formal) expansions with respect to qi as v and lu in (2.48) with new v'_ 
substituted for i/_„, 

We also assume that the new potential is real, i.e. U'* — U' and, correspondingly, 
L'* = L' . Moreover, we assume that there exists the resolvent M' of the new ope rator 
L' , that is its inverse operator. So we have also the second pair of equations in (2.30) 

MV = v'M , lu'M' = M uo'. (3.14) 

In order to solve the 9-equations for the new dressing operators in terms of the 
old ones it is convenient to consider, first, the 9-equations 

B^v'Wlo) = v'&W)^, ^(i/W-V) = v^W- 1 )^, (3.15) 
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that can be obtained from (2.51 ), (fP| ) and (3.5) noting that W and W~ x cancel 
exactly the additional S terms in ( |3. 5[ ) . In order to integrate explicitly these equations 
we write first that by (3.2) 

N N 

d^- 1 = S ( D i - a j) ( 3 - 16 ) 

3=1 



d 1 W = mY, c i 5 ( D i- h i)i 

3=1 

where 



N 



N 



i=i 



i=i 



N 



N 



(3.17) 



no,- 

i=i 



n ^ - 



z=i 



and ' means that the term I = j is omitted. Thus Q3.15 ) takes the form 

N 



3=1 



(3.18) 
(3.19) 



3=1 

where new operators were introduced whose kernels in the (p,q)-representation are 
independent on q and are, precisely, values of the dressing operators at some points 

v' b {p\ q) = v\p\ ib), uj' a (p; q) = uj'{p; ia - Pl ), (3.20) 

" a (Kq) = I/ (j ) ;i«)> w b (p;q) = u(p;ib- p t ). (3.21) 



Thanks to the composition law ( |2.3j ) and (3.7) it is easy to see that v a 5(D 1 — a) = 
v8{D 1 — a), S(D 1 — b)uj b = 5(D 1 — b)u>, etc. Just these equalities were used in (3.18) 
and ( 3.19| ) where in the r.h.s.'s thanks to them only ^-functions have kernels that in 
(p,q)-representation depend on q x . This means that we can rewrite (3.18) and (3.19) 

as 



N 



3=1 1 3 



0, 



TV 



3=1 



(3.22) 



(3.23) 



Since the kernels of v' , v, u>', u> and W in the (p,q)-representation tend to S(p) when 
qi — > oo, the expressions in parenthesis tend to / in the same limit and the equations 
can be explicitly integrated. Taking into account (2.28) we have 



JY 



3=1 



Di 



N 



/ = Wu>\I + Y i v« i 



3=1 



Dl -a 3 



(3.24) 
(3.25) 
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3. 2. Scalar products of the dressing operators 

Representations ( 3.24 ) and ( |3.25 ) obtained in the previous subsection are still 
undetermined, since they include the unknown multiplication operators v' b . and Lu' a _. 
In order to define them we need, first, to evaluate the scalar product of the dressing 
operators 

S = uj'v'. (3.26) 



Directly from this definition and thanks to (3.5) we have 



N 



JV 

E« 



where we used notation (3.8) for the residuum at a pole. In getting, for instance, the 
second term in the r.h.s. we used the relation lo'v'^D^ — a) = u)'{v'{D 1 —a))5{D 1 —a) = 
(S(D 1 —a))6(D 1 —a) = S a S(D 1 —a), which is obtained by noting that uj' can be inserted 
into the bracket since it has no pole in a. 

Thanks to ( 3.1 2| ) and (3.14) S must commute with Lq and Mq. The condition of 
commutativity of S with Mq formally is a consequence of the commutativity with 
Lq as their are inverse one to another. However, as we already mentioned, the 
composition law (|]^) is not necessarly associative and therefore both conditions must 
be imposed. Let us consider [Lq,S] — in the (p,q)-representation. We obtain 
from (2.3) and fl2.15| ) the equation [ip 2 — Pi(p x + 2q x )] S(p; q) = 0. Then S(p;q) 
must be a generalized function concentrated at p = and at p + 2^(qi) = and, 
consequently, a linear combination with arbitrary coefficients depending on qi of 
5(p), S (p + 2f sR-(qi)) and their derivatives up to a finite order. The derivative to 



be considered 



,( 



d 
9p t 



2i(Pi + Qi 



1-2- 



since it must annihilate [ip 2 ~ PiiPx + 2q 1 )]. 
This derivative when applied to 8 (p + 2£j}(qi)) can be equivalently written as g=- 
and, therefore, we conclude that the most general S commuting with Lq has a kernel 
which is a finite linear combination (with coefficients depending on qjj of the following 
distributions 



d 



-2i(p 1 + q 1 ) 



d 



gn 

S(p), Z n (p;q} = — S(p 



24i(qi)), (3.28) 



dp x dp 2 _ 

where n = 0, 1, Let us consider now the condition [Mo, S] = 0. One can check 

that only Yq and Zq commute with Mq, while Y n and Z n for n = 1,2,... do not 
commute with it and, precisely, the most singular terms in [Mo,l^] and in [MQ,Z n ] 
are, respectively, proportional to g- n -i S (iq2 — q?) S(p) and 5 (iq2 — qf) 5(p + 

2^si(qi)), which, of course, does not contradict the statement that Y n and Z n commute 
with Lq. Therefore S can only be a linear combination of Yq and Zq. But the term with 
Zq substituted in the l.h.s. of (3.27) produces a term Z\ that cannot be compensated 
by a term in the r.h.s. 

Thus finally we conclude that kernel of the operator S in the (p,q)-representation 
must be of the form 

S(p;q) = s( qi )S(p). (3.29) 

Now turning back to the equation ( |3.27 ) we see that in the (p,q)-representation the 
term in the l.h.s. as well as the second and third terms in the r.h.s. are proportional to 
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S(p) while the first term in the r.h.s. is proportional to S(p + 2£^(qi)) due to (2.52). 
So this term has to be equal to zero: 

[S,R'](p;q) ee ( S (-q 1 )- S (q 1 ))r'(q 1 )% + 2^(q 1 ))=0 



and ( 3.27 ) is reduced to 

N N 
3=1 3=1 



Taking into account that thanks to (3.6) and (3.26) at qi-infmity S(p;q) — > 5{p) we 
get 



or by ( |3.29D 



N 



»(qi) = i + E 



iS'„ 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



Finally, it is clear that for generic r(qi) (and, thus, r'(qi)) equation ( 3.30| ) gives 
that S a . = Sb =0 for all j and thus we proved that 

S = wV = /. (3.34) 

In other words we proved that like in th e case of decaying potential the dressing 
operators obey the first equality in ( [?.2£ ). It is necessary to mention that, on the 
contrary, as shown below, the second equality for the product v'uj' changes essentially. 



3. 3. Construction of Jost solutions and potential 

In order to proceed with the construction of the dressing operators v' and to' it 
is convenient to work in the x-representation and to introduce the corr espon ding 
s oluti ons <&'(x,k) and ^'(x,k) by means of the ' analog of the relations ( 2.31 ) and 
(2.32). Let us first consider the object f h = (Di — b)^ 1 ^^ appearing in the expression 
for v' in ( [3~24l ). By (|T|), ( |2~32[ ) and j23l| ) we have 



f h ( Xj x '. q ) = gg) J da e iWa+i 9l )-/(tt))«+ta(*i-« 1 ) >F ( a . j i6) a + iq j (3.35) 

where the so called Cauchy-Baker-Akhiezer function introduced in [p3 



■11 

^(x,k,k')= J dy 1 *( 2/ ,k)$(y,k') 



(k a -k' 3 )c 
V2 = x 2 



had appeared. Then by using (2.34) we obtain 



N , / N 

$/ ( x ' k )=IIk3 i - L *(^,k)+^c J $'(x,i6 J )^(x,i^.,k) 
i=i 1<2/ \ 3=1 



Analogously from lu' in (3.25) we have 



N 



k — ia 



JY 



*'(x, k) = I] j—t/ k ) - E k > ia 3) 



z=i 



3=1 



(3.36) 



(3.37) 



(3.38) 
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The function ^{x, k, k') obeys the 9-equations 
d 



^{x, k, k') - iTr^k - k') + T{x, k, -k')r(k'), 
k, k') = -inSCk - k') - T(x, -k, k')r(-k) 

ok 



(3.39) 
(3.40) 



= (see 
must be 



that are often considered as a generalization of equations (2.56). 

If (f>'(x,k) and ^'(a;,k) have no pole singularities, i.e. <&^A%) = ^ >' b . (x ) 
( |3.11 )) for all j, then the denominators in the r.h.s.'s of (3.37) and ( 3.38| ) 
compensated by zeroes of the expressions in parentheses. These 2N equations uniquely 
determine the 2N functions $'(x, ib.) and ^'(cc, ia^) and then the Jost solutions 
$>'{x, k) and ^>'(x, k) themselves. In this case it is easy to check that the corresponding 
dressing operators obey ( [3.12 ) that proves the statement given there for the absence 
of poles. Turning back to the generic situation, in order to complete the construction 
of <&'(£, k) and \I>'(a;,k) we evaluate the values &'(x,ibj) and ^'(x, ia^) appearing 
in (3.37) and (3.38) by using the equality lo'u' = I derived in the previous section. 
Once constructed $'(x,\i) and ^'(a^k) in terms of <I>(x,k) and ^(a^k), by using the 
asym ptoti c limits lim-j^^oo k) = lirn^^oo £(x, k) = 1, easily derived from (2.40) 
and ( 2.41 ), we must verify that %'(a;,k) and £'(:r,k) are polynomially bounded at 
space infinity and, then, that v' and lo' belong to the considered space of operators, 
which guarantees the correctness of the above derived equalities and, in particular, of 
( 3.34 ). We assume this behavior and below we shall prove it in some special cases. 

It was already mentio ned that the first equality in ( 2.28 ) can be rewritten in 
terms of $ and as ( 2.44 ). The same orth ogonality condition for $' and '5' follows 
from ( 3.34 ). Inserting in it ( 3.37 ) and ( [3.38] ) we get 



I N 

dx 1 I ^(x,'k + p) — y y^c j J r (x,'k + p,ia^'(x,ia^ / 

\ 3=1 



N 



x j $(x,k) +^c J &{x 1 ib J )T(x:ib ]1 1 k) ] =2tt5(p). 



(3.41) 



Thanks to the above assumption it is clear that the integral exists in the sense 
of distributions. Now we multiply this equality by *&(x', \i)&( x" , k + p) where 
x 2 = x 2 — x 2i integrate it with respect to ks^ and p and use ( f2~45| ). Then the 
r.h.s. becomes equal to Att 2 S(x' 1 — x'{). Correspondingly, the l.h.s. at a generic point 
k must have zero ^-derivative and the residua of poles at k = ioj and k = ibj must 
be equated to zero getting 



jv 



V'faiaj) = ^2T'{x,\a j ,\b l )c l ^!{x,ib l ), 



i=i 



N 



^'{x,ibj) = -^$(x,iaj)cj.F'(x,iaj,i&,-)> 
i=i 

where in analogy with ( 3.36[ ) we introduced 



F'(x,k,k') 



d Vl *'(y,k)$'(y,k'). 



(3.42) 
(3.43) 

(3.44) 
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A' 



^'(.T,ia J ,i6 ; )c ; *(x,i6 z )$( x, \a m )c ra T (x, ia m , \b k ). (3.45) 

1,171 = 1 



Now by (3.44) we rewrite this as 
d 

— F(x,i aj ,ib k ) 



A 



or as 



l,m=l 



8 .(JF')-l = A c ^ 5i 



(3.47) 
(3.48) 



dx 1 dx 1 

where we introduced matrices T and T' with elements 

Fjk = F(x, ibj,ia k ), T' jk = T'{x, ia 3 ,ib k 

and the two diagonal matrices c = diagjci, . . . c/v} and c = diagjci, . . . c/v}. Let C 
denote an a^-independent matrix. Thus we get 

T' = c- 1 {C + T)- 1 c- 1 (3.49) 

and inserting T' in ( 3.42 ) and ( 3.43| ), correspondingly, 

1 N 

*'(x,ia j ) = T J £(C + F)j l 1 *(x,ib l ), 



"o 1=1 
N 



(3.50) 
(3.51) 



°j i=i 

In order to avoid singularities we have to impose the condition that det(C + ^F(x)) 
has no zeroes on the x-plane. At the end of this Section we show that in the cases 
N = 1 and N = 2 this condition is formulated as a condition on admissible matrices 
C. For generic N we assume the existence of such matrices C that guarantee absence 
of zeroes of this determinant. In additi on on e has to check that for these C the matrix 
T' with elements T'- k obtained from (3.44) by taking k = ia,j and k' = ib k and by 
s ubst ituting, respectively, ( |3.5C| ) and (|3.5l| ) for &{x,ib k ) and fy'(x,iaj) coincides with 
( |3.49 ) . This leads to the equality 



(3.52) 



(aj-bk)c 

i.e. all elements of this matrix given by ( 3.49| ) have to obey conditions 

lim ?- k (x)=0. (3.53) 

It is easy to check these conditions in the cases N = 1 and N = 2, but for generic N 
we have to assume their validity. 
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Now from (3.37) and (3.38) we get for the residua of <&' and 

N N 

i=i i=i 

We can also insert the r.h.s.'s of ( 3.5CD and ( 3.5l| ) into (3.37) and (3.38) getting 



N 

m— 1 
N 



k - ib r 



k — ia 



N 



$(x,k)- E Hx,i aj )(C + ^)j^(x,ib h k) 



3,1=1 
N 



(3.55) 



(3.56) 



The above formulas can be rewritten in terms of ratio of determinants as it was done 
for the KPI-casc [Ell. They can be considered as a special case of the construction 
given in p2| for generic eigenfunctions, i.e. not necessarily Jost solutions. On the 
contrary, here we deal with objects that indeed can be called Jost solutions as their 
departure from analyticity is under control and follows from ( |3.5| ) by the ' analog of 
the relations (2.31) and (2.32) and equalities (3.10): 



N 



d k <Z>'(x, k) = -k)r'(k) + bT]T K 3 (s)*(k - iaj). 



j=i 



iV 



flk^'Ca;, k) = -*'(a;, -k)r'(-k) + i^E *t, (z)«(k - 



(3.57) 



(3.58) 



j=i 



where $^(a;) and %.{x) obey (flij). Thanks to (|3~TT|) and ( |3~55| ), ( |3~5^) we get 
explicitly 



JV 



0) = ^ E *( x M){c+?)tA 



N 



KM 



(3.59) 



(3.60) 



/,m— 1 



where we used notations ( |3.17| ) and (3.48). In particular, the above discussed case 
where $'(a;,k) and ^'(a^k) have no poles corresponds to matrix C = 0. These 
transformations do not add solitons and can be viewed as transformations of the 
continuous spectrum. 

Now we can derive an explicit formula for the potential v! with spectral data 
(|3.3[). Inserting in (3.55) the two leading terms of the asymptotic expansion (2.- 



and ( 2.48 )', i.e. x( x ,k) = 1— (2ik) 1 X-i( x ) + • • • and its analog for we get thanks 
to (^3lT 

N 

X-i(^)=X-i(^) + 2E( <I j- b j)- 29 x I logdet(C+JF(x)), (3.61) 

3=1 

so that for the potential we have 

u(x) = u{x) - 2d 2 Xl logdet(C* + T{x)). (3.62) 
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It is easy to check directly that if C'(x,d x ) denote the operator of the type 
(1.1) with potential u'(x) substituted for u(x) then ( [3.55] ) and (3.56) obey equations 



C(x, d x )$'(x, k) = and its dual if and only if 

d X2 C = 0. (3.63) 



Formula (3.62) is well known in the literature |2l], [2^], but nevertheless conditions 
that guarantee the regularity of the new potential are, to our knowledge, absent. 
This situation is completely different from the case of the Nonstationary Schrodinger 



Equation, where conditions of regularity were given in Eq|. If all aj and bj in (3.2) are 



real, it is enough to choose the matrix C real in order to get a real potential. However, 



in the generic situation when (3.4) is satisfied, the conditions to impose to C in order 
to have a real potential are unknown. Even more complicated are the requirements 
to impose to the matrix C in order to obtain regular potentials. 

We made a complete analysis in the cases N = 1 and N = 2 when the aj and bj 
are real. Then the matrix C must be real. In the case N = 1 the matrix C reduces 
to a real constant and the regularity condition is (ai — b\)C > 0. The case C = 
corresponds to a non solitonic situation and, in particular, when the original potential 
u(x) is equal to zero we get that u'(x) — also. In the case N = 2 a complete analysis 
is elementary but rather lengthy. Here, we give only the final result, more details are 
presented in the Appendix. It is convenient to introduce the following constants 

Dij = dj Ji+i,j+i, Jij = — , (3.64) 

CLj Oi 

where i and j are defined mod 2 and to choose, for definiteness, ai < a 2 and bi < 6 2 . 

In order to have a potential describing two solitons superimposed to a background 
at least two must be different from zero. 

Then the potential is regular and has four rays if and only if for i = 1, 2 

in the case dct J < 0, J\\J%i > 

£>u<0, D i)i+ i>Q, detC<0 (3.65) 
in the case det J < 0, Jn J22 < 

Da < 0, D M+1 > 0, det C < (3.66) 
in the case det J > 

Da > 0, A,i+i < 0. (3.67) 

or 

D u > 0, D iA+1 < 0. (3.68) 
The potential is regular and has three rays if and only if for i = 1 or 2 and j — 1 or 2 
in the case det J < 

D u = D jtj+ i = 0, A+m+i < 0, D j+hj > (3.69) 
in the case det J > 

D u = D j>j+1 = 0, A+m+i > 0, D j+1J < 0. (3.70) 
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If three are zero the potential contains only one soliton. 

Notice that the two soliton solution u'(x) can be considered as a surface in the 
(x, u') space depending on 8 parameters, i.e. — 1, 2), <Zj and 6, (i — 1, 2). Since 

in the space of parameters the surfaces Dij — and det C — separate singular from 
regular solutions the soliton solution is not differentiable on these surfaces, which are 
therefore bifurcation surfaces according to the usual definition in catastrophe theory. 
Therefore, we expect that the solution, as a geometrical object, would be structurally 
instable at these values of the parameters. 

The possible behaviours of the potential at large distances in the x-plane, 
according to the different choices of the parameters, are richer than in the KPI case. 
For instance, in the case det J > 0, if we call x\ + hx2 — const, the direction of a 
soliton, the directions at X2 — ±oo of the two solitons when the matrix C is full are 
given by h = ai + a2 and h = b\ + &2> when the matrix C is diagonal by h = a\ + b\ 
and h = a 2 + b 2 and when the matrix C is off-diagonal by h = a\ + b<i and h = a 2 + b±. 
In addition if only one element of the matrix C is zero the four rays of the solitons are 
directed along four different directions. These different behaviours are in agreement 
with the previous comment on the structural stability of the soliton solution. In the 
appendix we give a detailed description of the potentials in all regular cases. 



4. Completeness relation and resolvent 



We proved in section 3.2 that, like in the case of decaying potentials, the scalar 
product (3.34) wV is equal to I or, in other words, that v' is right inverse of oj' . 
On the contrary we expect, due to the presence of discrete data in the spectrum, that 
the second equality in ( 2.28| ), the so called completeness relation, is modified by an 
additional operator P' as follows 

v'u' + P' = I. (4.1) 



From (3.34) we get directly that 



PV = 0, 



JP 1 = 0, 



p' 



p' 



(4.2) 



i.e. P' is an or thogo nal projector. By substituting in v'u)' the values of v' and oj' given 
in ( 3.24 ) and ( 3.25| ) we get an expression of P' in terms of the residua of v' and oj' 



N 



D 



i , sr / 1 



3 = 1 



D, 



(4.3) 



that in the x-representation can be written explicitly as 
P'(x, x'-q) = - sgn^ - x[)S(x 2 - 4)e-^-^ 

JV 

x £Wfe - b m ){x x - x'^'ix.ib^.ix') 

3=1 

+9((q l -a M )(x 1 -x , l ))& a (x)V , (x , ,ia j )}. 



(4.4) 



This equality can be simplified by means of (3.54) 
P'{x, x';q) = -S(x 2 - x' 2 )e- qi{xi - x '^ 



N 

x ^{x 1 ib J )c J C ]l c^'{x\ia l ){9{q 1 - b M ) - 6{q x - a m )}. (4.5) 

3,1=1 
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Then we see that P' is equal to zero outside the largest interval in qi with extremes 
a m , (j, I = 1, . . . , N) and that dP' (see (2.6)) is equal to zero if qi is different from 

and bjsft. Again it is necessary to emphasize that all these results are obtained on 
the basis of the assumption of polynomially boundedness of the x' ancl £' formulated 
in the previous section. 

Also the dressing formula for the resolvent M', inverse of L', must be corrected 
with respect to the bilinear representation (2.29) adding an operator ml as follows 

M' = v'Mqlj' + ml. (4.6) 



Thanks to ( 3.12 ) we have L'M' = L'u'Mqlo' + L'ml = v'ui' + L'ml, where we u sed the 
fact that LqMq = I. Thus in order to obey L'M' = I the operator ml by (4.1) has to 
obey the equality L'ml = P' . 

Let us consider here the case N = 1 and let a = a±, b = b\ be real and let us 
label the corresponding quantities with 1. We have then 

M l =v x M Q u l +m l (4.7) 



and, in this case, it is easy to derive from (4.3) that 

1 



' (D 1 - a)(D 2 -(a + b)D 1 + ab) 

1 



1 ' b {D 1 -b)(D 2 -(a + b)D 1 +ab) ^ 



(4.8) 



In writing the kernel of thi s operator in the ^-representation we can use ( |3.54|) that 
in this case thanks to ( 3.17 ) takes the form 

$i,a(z) = ^(x/ib) ¥ lf6 (a:) = -cV^ia), (4.9) 

where 

c={a-b) 2 C ll . (4.10) 
We get a formula similar to ( |4.5| ) , that is 

mi (x, x'-q) = (6( qi - a) - 9( qi - b)) [9{x 2 - x' 2 ) - 9(~q 2 + (a + b)q ± - ab)] 

xce-t^-^^ix, ib)^ x (x\ ia). (4.11) 

From this equation we have L\m,\ — P\ as required. Let us, however, stress that 
we would obtain this equality for any other x-independcnt second term in the square 
bracket and that the above specific choice is necessary in order to avoid an exponential 
growth of mi(x,x';q) at space infinities, as easily follows from relations (3.50) and 
(Hi]). For generic N the construction of ml and the study of its asymptotic properties 
are essentially more cumbersome. 

Let us discuss here some properties of the resolvent M\ constructed for the case 



N = 1 as they follow from (4.7) and ( ]4.q )-(4.11). We see that, in comparison with 
the resolvent M ( |2.29| ) of a decaying potential, Mi has an additional discontinuity at 
92 = (a + b)qi — ab due to the term mi(i, x'; q). This discontinuity is not compensated 
by the first term in ( |4.7| ) and its presence in the resolvent Mi is a characteristic 
manifestation of the solitonic content of the potential u\(x). The term mi(x,x';q) is 
zero when q\ is outside the interval (a, b) and therefore discontinuous also along the 
lines qi — a and qi = b on the g-plane. These discontinuities can be compensated by 
the pole behavior of the dressing operators v\ and u>\ in the dressed term viM q uj 1 and 
need a detailed study. 
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From (3.5) using the Cauchy-Green formula and the asymptotic behaviour of V\ 
and u>\ we have 



, + I/^ 1(Rl „ l) <., ) + _L_„ v , 



where notation (2.5) for the shift was used. Then inserting them into V\MqU)i we get 



v 1 M uj 1 = Mq 



I f iW[fili),^] u W + M 1)discr , (4.14) 
~ J z i 



where 



M i,discr = "i,»n w i,a + "tin E^i.b- ( 4 - 15 ) 

In the x- representation we get for the contribution of the discrete part of the spectrum 
to vxMqUx, thanks to (JyJ|), fl2~2(i| ), ( |3~To| ), and ( ^54| ), 

M 1AiscI (x, x'; q) = c^(x, ib)* x {x', ia)e-^-') [T a (x - x'; q) - T b (x - x'; q)], (4.16) 
where 

i£(ia)x p — U(a+iqi)x 

r a (^; <z) = -5— /da[% 2 - g 2 - a 2 ) - 0(a; 2 )] -. (4.17) 

2ir J q 1 — a — la 

The singularities of r a can be explicitly extracted getting 

r a 0; l) = ^reg^! l) + V u, S in R { x \ ?)) ( 4 - 18 ) 

where 

p , x %2) J, f xi + 2ao; 2 
r a)reg (x; g ) = ^-|l-erf^-^- 



•\AF-92 

/ da* — r - (4.19) 

27ri J a + i(q 1 — a) 



and 



T a sing (:r; q) = -9fa)9{q x - a) + ^L^) arctan g 2 (4 2Q) 

' & 7r q x — a 

We conclude that, for generic g 2 , the discontinuities along the lines q± = a and q\ = b 
of Afi.discr(a;, x'; q) cancel exactly the discontinuities along the same lines of mi, but 
e -i x Y a<S i ng (x;q) and e~ qx Tb, s ing(x;q) in the neighborhood, respectively, of the points 
(a, a 2 ) and (6, b 2 ) in the g-plane are ill defined due to the arctan. 

One of the main advantages in using the resolvent is that the Green's functions can 
be obtained as specific reduction with respect to the parameter q of the resolvent. In 
particular, if we are interested in considering a smooth perturbation it 2 (a:) (decaying 
at space infinity) of the potential U\{x) constructed above, then the Jost solution 
$(x,k) of the perturbed potential u = ui + it 2 can be obtained as a perturbation of 
the Jost solution $x(x,k) by means of the following integral equation 

3>(x,k) = $ x (x,k) +Jdx'G 1 (x,x',k)u 2 (x')$(x',k), (4.21) 
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where (cf. ( |2~38j )) Gi{x,x',k) = e^-^M^x' ; g) | ^ _^ ^ . Taking into account that 

thanks to (2.23) £$(k) = (kg, — k^+k^) we get that the above mentioned discontinuity 

coming from mi becomes a discontinuity across the hyperbole (ka — 9 ^) 2 — kjjj = 

(^tO 2 wn i cn ues m the k-plane outside the strip a < < b if a < b or & < k$> < a 
if 6 < a. On the other side the term in G\ coming from mi is equal to zero outside 
this strip. Therefore, according to our previous discussion, only the discontinuities 
at k = ia and at k — ib are left. Thus in the case of the heat equation (at least for 
N = 1) the Green's function of the Jost solution has no an additional cuts in contrast 
with the case of the nonstationary Schrodinger equation [TBI , but, anyway, due to 
the special singularities at k = ia and at k = ib, the definition of the spectral data also 
for the case of a perturbed one soliton potential is not standard and needs a detailed 
analysis. This problem will be faced in a following work. 
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Appendix 

The potential u'(x) describing two solitons superimposed to a generic background has 
at large distances in the x-plane a solitonic one dimensional behaviour along some 
rays. Any specific such ray can be represented by an equation xi + hxi = const, with 
a given h and by specifying if the rays is pointing to X2 = +oo or to 12 = —00. We 
choose, for definitencss, 02 > a\ and 62 > b\. 

Let us, first, list the cases in which the potentials has four rays. 

It is convenient to rename the parameters at and bt (i = 1, 2) as ctj (j = 1, . . . , 4) 
in such a way that 

a\ < a 2 < a 3 < a^. (A.l) 

If Cij and det C^Owc have the following rays 



h 


versus 




±00 


a 2 + Q!4 


±00 



(A.2) 



Note that the rays are two by two parallel. 

According to the remark we made at the end of Section 3 we expect that when 
an element of the matrix C or det C is taken to be zero the asymptotic behaviour at 
large distances will change discontinuously. If only one element of the matrix C or 
det C is zero we have four rays at large distances but they are no longer two by two 
parallel. If two elements not belonging to the same diagonal of C are zero we are left 
with only three rays. If three elements are zero the solution reduces to the one soliton 
solution and if the matrix C is zero the solution does not contain solitons. 

The details are the following. 

If, in the space of parameters of the soliton solution, we move to a bifurcation 
point by taking Cu = for i = 1 and/or 2 (det C ^ 0) the four rays are obtained 
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by transforming discontinuously among the rays listed in the previous table (A. 2) the 
following rays for i = 1 and/or 2 



into the rays 



h 


versus 


a\ + a 3 


R' 1+i J n oo 


OL2 + OL4 


{-) l+L J u oc 




h 


versus 


Oil + OL4 


R i+i J n oo 


®-2 + CC3 


J n oo 



(A.3) 



(A.4) 



If Oj,i+i = for i — 1 and/or 2 and det C ^ the four rays are obtained by substituting 
for i — 1 and/or 2 the following rays 



h 


versus 


ol\ + a 3 


(-Too 


OL2 + OL4 


(-) 00 



among the rays listed in the table ( A. 2 ) with the rays 



h 


versus 




(-)'oo 


OLZ + ot A 


(-) oo 



(A.5) 



(A.6) 



If det C — 0, det J < and Jn J22 > and all dj the four rays are 



h 


versus 


OL\ + OL2 


— Jnoo 


OL\ + Qf3 


+ JllOO 


Oi2 + OLA 


— Jnoo 


OLz + 014 


+ JllOO 



(A.7) 



If det C = 0, det J < and J n J 2 2 < and all CV,- ^ the four rays are 



h 


versus 


ai + as 


— Jnoo 


Oi\ + 0.4 


+J1100 


012 + Oi 3 


+ JllOO 


OL2 + OL4 


— Jnoo 



(A., 



Let us, now, list the cases in which the potential has three rays. 
If Cu = and C^i+i = for i = 1 or 2 the three rays are 



h 


versus 


ai + a 2 


+ Ji+l, i+l Jj+l/jOO 


at + b l+ i 


J i+ i,i + i<X> 


a-i+i + &i+i 


(~y Ji+i,ioo 



(A.9) 



If Cu = and Cj+i j = for i = 1 or 2 the three rays are 



h 


versus 


h +b 2 


— Jj+l,i+l Ji,i+lOO 


a i+ i + bi 


Ji+i^+iOO 


o-i+i + 6-i+i 


(~Y Jia+lOO 



(A.10) 



Finally, let us note that the corrections to this one dimensional solitonic behaviour 
are exponentially decaying, in contrast with the KPI case, where the corrections are 
rationally decaying at least in some regions of the plane |Q, 
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